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Abstract 

This paper is devoted to the estimation of a vector parametrizing an energy 
function associated to some "Nearest-Neighbours" Gibbs point process, via the 
pseudo-hkeHhood method. We present some convergence results concerning this 
estimator, that is strong consistency and asymptotic normahty, when only a 
single realization is observed. Sufficient conditions are expressed in terms of the 
local energy function and are verified on some examples. 



1 Introduction 



Gibbs point processes first appeared in tlie tlie ory of statis t ical p hysics. Historical aspects 
of the mathematical theory are covered briefly in iKallenberd (|l983) . The importance of the 
Gibbs point process as a model building principle became widely recognized through these 
works. Indeed, the class of Gibbs point processes is interesting because it allows to introduce 
and study interactions between points through the modell ing of a.n asso ci ated potential fu nction. 
This resulting gain explains their use in statistical phvsics lR.uellel ()1969l ) . iFevnmanI ()1972l ) (when 
taking interactions between molecules in models of dilute gases into account) or in ecology (when 
analysing competitions between plants). Within the mechanics statistics framework, Gibbs states 
are defined as solutions of the well known e q uilibrium equations refered to Dobrushin-Lanford- 
Ruelle (D.L.R.) equations iDobrushinl (|l969l ^. iLanford and Ruelld (|l969l ). One way to introduce 
Gibbs point processes consists in using a family of local specifications with respect to a weight 
process. The Preston's theorems (jPrestonl (ll97y)) used precisely this approach in order to give 
sufficient conditions on local specifications for the existence of Gibbs states. 

Many proposals tried to estimate the potential function from the available point pattern data 
generated by some Gibbs point processes. If the potential belongs to a parametric family model, 
the most well-known methodology is the use of the likelihood function. The main drawback of 
this approach is that the likelihood function contains an unknown scaling factor whose value 
depends on the parameters and which is difficult to calculate. The first class of models on which 
the estimation of th e maximum likelihood has b een undertaken is the class of pairwise interac- 



tion point processes. iQgata and Tanemura (|l9M) 

developed th e maximum likelih ood estimation 

method based on numerical approximations of the likelihood. iPenttinenI ((198J) used a similar 
approach while applying a Monte Carl o method in a way to solve th e likelihood equation by the 
stochastic Newton- Raphson algorithm. iMoveed and BaddelevI (|l99ll ) proposed another iterative 
procedure for estimatin g the maximum likelihood est i mator. For m aximum likelihood by Markov 
chain Monte Carlo, see lGever and Thompson ( 19921 ) . iGeven ([19991) and for U. L.A.N, conditions 



for maximum likelihood estimator, see 



Mase 



Ji 19921). 



An alternative approach consists in avoiding 



to optimize the likelihood function (because of the scaling factor problem) and introducing a 
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pseudo-likelih ood function i nstea d. This idea originated from iBesad (|1974f ) in the study of lat- 



tice processes. 



rile 



process, w 
cesses, see 



Besag et al. J 1982 ) further considered this method for pairwise interaction point 



Jensen and Mollei { 1991 ) generalized it to the ge n eral cla, s s of G ibbs point pro- 



Masd (|199a ). lMasel (119991 ) . iJensen and Kiinschl (|l994r ). iGuvon 1991) for asymptot ic 
properties. A third way is the Takacs-Fiksel estimation method (|Takacsl (|1986l ) iFiksell (|l988l )). 
which relies on a characteristic property of Gibbs processes using Palm measur e. Asymptotics 
properties of Takacs-Fiksel estimator are studied in lHeinrichI (|1992l ). lBilliod ( 1997). A compariso n 



of these different procedures applied to the Strauss n iodel is presented in 



Diggle et al 



(1994) 



The non parame t ric se tting has been undertaken b v iGlotzl and Rauschenschwandtneii ()l98ll ) 



and 



Diggle et al 



( 19871 ) (and the references therein) . iHeikkinen and PenttinenI (|l999l ) proposed 
a semiparametric estimator based on Bayesian smoothing techniques. A general review of the 
pr oblem of statistical inferenc e on sp atial point processes can be found in the recent monograph 
of 



Mgiller and Waggeoetersenl (|200 



The present study is devo ted to "Nearest-Nei ghbour" Gibbs point models by c ombining 



stocha s tic geometry arguments (IStovan et al.l (11993) 



1988 ) , lEdelsbrunneii (|l988l ) , iBoissonnat and Yvined (|l99a ) ) . Such models are introduced by 



and c omputational geometry ones (IPreparata and Shamos 



Baddelev and Mgilk 



19891 ) where the neighbourhood relation depends on the realization of the process. Sufficient 



conditions ( expressed in terms o f the energy function) fo r the existence of such processes are 



proposed in 



Bertin et al 



(1999b|) and 



Bertin et al 



( 1999al ). where some examples are also pro- 



posed. The main one is a pairwise interaction point process where the neighbourhood relation 
corresponds to the (slightly modified) Delaunay graph of the realization of the process. 

In this paper, we study a pseudo-likelihood estimator for such processes. More precisely, our 
framework is restricted to stationary Gibbs point processes based on energy function related 
to some graph (for instance the Delaunay graph) such that the energy function is invariant by 
translation and such that the local energy function is stable and quasi-local (or local). The main 
results of this paper are convergence results (strong consistency and asymptotic normality) of 
maximum pseudo-likelihood estimators in this framework. These results are obtained when only 
a single realization is observed. Sufficient conditions are expressed in terms of the local energy 
function (which makes the results quite general) for some large family of parametrized energy 
functions. Among the different parametrizations, the exponential family is considered. 

The paper is organized as follows. Section |21 is devoted to some background on Gibbs point 
processes and to the description of our framework. The statistical model and the pseudo- 
likelihood method are presented in Section |21 Consistency and asymptotic normality of the 
maximum pseudo-likelihood estimator are respectively proved in Section |1] and Section [SI Fi- 
nally, the different sufficient conditions ensuring convergence results are verified on some ex- 
amples in Section |21 A short simulation is presented to check the effectiveness of maximum 
pseudo-likelihood estimator. 
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2 Background on Gibbs point processes 

2.1 Gibbs point processes 

We define B, to be respectively the Borcl cr-field and the bounded Borel boolean ring. 
Let Q. denotes the class of locally finite subsets of IR'^. In particular, an element 99 of fi, 
also called configuration (of points), could be represented as 99 = 'YlieYN^^i which is a simple 
counting Radon measure in IR*^ (i.e. all the points Xi of IR*^ are distinct) where for every A G 
B , 5x{^) = 1a (3;) is the Dirac measure and Ia(-) is the indicator function of a set A. This space 
is equipped with the vague topology, that is to say the weak topology for Radon measures 
with respect to the set of continuous functions vanishing outside a compact set. We also define 
the fj-field J- spanned by the maps ip — > > ^ ^ where V3(A) corresponds to the number 
of points of 93 in A due to the Radon measure representation of (p. The set of all configurations 
in a measurable set A C IR*^ will be denoted by Oa and the corresponding a-field J-\ is similarly 
defined. Furthermore, for any A G Bi,, 

= (j^a,^a) X (J^Ac,:^Ac) 

where A'^ = IR'' \ A denotes the complementary of A in IR^. Finally, ilf denotes the class of all 
finite subsets of JR'^. 

A point process on IR'' is a valued random variable, denoted by with probability distri- 
bution P on {Q,J^). and the intensity measure Ap of P is defined as a measure on B such that 
for any D E B 

Ap{D) = [ ^{D) P{d^). 

In the stationary case, Ap{D) = \pv{D) where the constant \p is called the intensity of P and 
V is the Lebcsguc measure on IR'^. 

A Gibbs point process is usually defined using a family of local specifications with respect to 
a weight process (often a stationary Poisson process with distribution Q and intensity Ag = 1). 
Let A be a bounded region in IR'^. For such a process, given some configuration 99 a^ the 
conditional probability on A is of the form, for any Y E T : 

JlK{ip, Y) = [z~(^J^ ^""P ^"'^ W<PAc)) ly (V' U ¥'ac)(3a(#)} IrM, 

where 

Zk{^)= i exp(-y (V'|(yCAc))(5A(c^V') 
7 Ha 

is called the partition function and R\ = G J7 : < Z^{}p) < 00}. 

Whereas the finite energy function V (ip) measures the cost of any configuration, the local 
energy V {tp\(p) is defined as the energy required to add the points oi tp in (p : 
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Let us notice that when ^p reduces to one point x, we denote by a shght abuse V {x\ip) instead 
of V {{x}\ip). It is weh known that the collection of probability kernels (nA)AeBb satisfies the set 
of co mpatibility and measurability conditions which define a local specification in the Preston's 
sense (|Prestonl (jlOTg)). The main condition is the consistency : 



A' 



Ha for A C A. 



Notice that some conditions are needed to ensure the existence of a probability measure P 
with respect to any local energy V and any weight process that satisfies the so-called Dobrushin- 
Lanford-Ruelle (D.L.R.) equations : 

P(y|:rAc)(v?) = nA((^, Y) for P a.e. ipen for any A G and F G J^. 



For the general theory of G i bbs point processes , the reader may refer to iKallenberd (|l983r ) 
Daley and Vere-Joned (|l988l ): 



Stovan et alJ (|l995|) and the references therein. 



2.2 Campbell and Palm measures and Glotz Theorem 

The reduced Campbell measure Cp ofP is a measure on such that for any D £ B and 

any Y ^ T 



CUd X Y) 



in J D 



When some measurable function h from IR'* X Q on is given, the following equation is often 
called the refined Campbell theorem 

/ ^h{x,ip - 6^)P{dip) = / h{x,ip)Cl{d{x,ip)). 

If the intensity measure Ap is a— finite, then for Ap- a. a. x G IR*^, the distribution P'^ on {^l,J-) 
exists. It is unique for Ap- a.a. x G IR*^ and such that 

CUd xY) = [ pI(Y) Ap{dx) for any D£B,Y e T. 

Then PI. is called the reduced Palm distribution of the point process P with respect to point x. 
Intuitively, the Palm distribution Px is the conditional probability of configurations of the point 
process given that the point x belongs to the realization (/?. Therefore, we have 



m.\n 



h{x,ip)Pxid(p)Ap{dx). 



When the process is stationary, one may apply the previous equation by rep lacing the intensity 
measure Ap by its expression in this case Xpi', and in this framework, Glotzll |l98Q) proved that 
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P & Go if and only if the reduced Campbell measure Cp is absolutely continuous with respect 
to X P and : 

dCp , . , dP- 



{x, y?) = Ap— f ((^) = exp {-V {x\^)) 



d{vxP)^ dP 

where Xp =J exp {—V {x\ip)) P {dip) is the intensity of the process P. In the particular case 
n 

when V {x\ip) = 0, the point process corresponds to the stationary Poisson process Q. We know 
from the Slivnyak's theorem that Q'^ = Q which is one way of characterizing such process. 

2.3 Description of some Gibbs models 

This paper is mainly dev oted to the statistical s t udy o f some nearest-neighbours Gibbs point 



processes first introduced in iBaddelev and Molleii ((19891). More precisely, we are interested in 
models based on energy function of the form 

^M = E E ^^'H^,^), (1) 

k=l gG-Deife((/3) 

where Delkif) is the set of clique of order k of the Delaunay graph defined just below. For some 
(/? € $1 in general position, one defines Del^{ip) by the unique decomposition into triangles ij) 
in which the convex hull of the circle C{il)) does not contain any point of </9 \ ■0. The Delaunay 
graph is then defined by the set of edges : 



In order to ensure the existence of such Gibbs state in IR'^. lBertin et al.l ()l999bl ^ prove that the 
local stability and quasilocality properties (only expressed in terms of the energy function) are 
sufficient conditions of Preston's Theorem. 

Without any additional modification, the previous model does not satisfy the previous as- 
sumptions. We then introduce some subgraphs. First let us denote, for some triangle ■0, by D{iIj) 
the diameter of the circle circumscribed of ^ and by /? (ip) the smallest angle of ip. 

Definition 1 Given any Pq G]0,7r/3], we introduce the following particular subset of Del^ (ip) : 

De4°p iv) = {V G Dek (v) ■ P W > f3o} • 
The (3-Delaunay graph of order /?q of any configuration ip is the Delaunay subgraph defined by : 

The model obtained by replacing the original Delaunay graph by the /?-Delaunay subgraph of 
order /3q in satisfies the previous sufficient conditions of Preston's Theorem. From now on, 
this model is called the /3-Delaunay model. In this spirit, some other models may be defined (see 
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e.g. iBertin et alJ (|l999b|), iBertin et alJ (ll999a|)) but we advice the reader to keep in mind the 
/3-Delaunay model as the main example in order to illustrate the statistical results developped 
in this work. 

The framework of this paper is restricted to stationary Gibbs point processes based on energy 
function related to some graph, denoted G2{f) for some finite configuration ip (Gkif) representing 
the set of cliques of order k) , of the form 



^max 



k=l lceGfc(</p) 



(2) 



J^max 



E when^xW^^W 

k=2 IcGGfcH 



and satisfying Assumptions Ei, £3°^^ or more generally £2^°^, E3 defined by : 



El V (•) is invariant by translation. 

E2°'^ Locality of the local energy : there exists some fixed range denoted by D such that for 
any 99 G one has 

v{<d\<f) = y (o|v9ni3(o,L>)) . 

^2^''^ Quasi-locality of the local energy : there exists a nonnegative function e vanishing 
asymptotically such that for any (/? G O one has 

\V (Olv?) - V (0|v9 n 0(0, D)) I < e{D). 

E3 Stability of the local energy : there exists K >{) such that for any € 0, 

V{{)\ip) > -K. 

This framework includes some classical point processes such as : 

- models based on the usual complete graph G2if) = ^2(9')) with pairwise interaction 
function satisfying a hard-core or inhibition condition and with finite range. 

- /c-nearest neighbours models w ith pairwise interaction function bounded and with finite 



range (see 



Bertin et al 



(1999d)). 



Widom-Rowlinson or area interaction model. 



3 Statistical model and inference method 

3.1 Statistical model 

We consider Gibbs point processes with energy function V {■■,6) parametrized as follows 
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As a statistical model we consider a parametrized version of ((21) where the different u^^\^, (f) 
depend on a vector parameters O^''^ and then denoted from now by u^^\^; if, 0^^^). It is assumed 
that the vector of parameters = {9i, ... , Op+i) = (^9^^\6^'^\. . . , e where is an 

open bounded set of IR^^^. 

Our data consist in the realization of a point process with energy function V (•; 9*) in a do- 
main A C IR'^ satisfying Assumptions Ei to E3. Thus, 0* is the true parameter to be estimated. 
The Gibbs measure will be denoted by P0* . From 0, we obtain easily the energy to insert a 
point X in a configuration ip. 



Krr 



V{x\ip;e)= J2 < 



k=l 



UeGfc(<^uM)\Gfc(¥^) i&Gk{'p)\Gk{^\J{x}) 



(3) 

Theoretical results presented in the next sections are valid for a general energy function 
V {-^O). But among this class of models, we will focus on energy functions described by ^ and 
such that 

The energy can be rewritten 

k=l i^Gkiv) fc=l 
where for any finite configuration ip 

u^'\^)= E ^^"H^;^) and wM = (^i(^),...,nj,+i((^))= )M) 

For two finite configurations (p and ip, by denoting 

u^^\^\ip) = u^''\ijj U ip) - u^''\ip) and u{7p\ip) = (u^^'Xiplip), . . . ,u^^"'''-\ip\ip)J (5) 
we have for any point x 

V{ip;e)= Y ^^^^ '^^^H^) = ^u{ip) and F(x|(^;6>)= ^ ^ '^^''^ i^l^P) = ^ ^(^\^) ^ 

k=l k=l 

(6) 

where u{x\(p) = {ui{x\(p , . . . , Up-^i{x\(p) = {u^^\x\(p) , . . . , u^™-"-"" {x\(p)) . The local specification of 
the Gibbs point process associated to an energy function defined by © belongs to an exponential 
family. 
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3.2 Pseudo-likelihood 



As precised in the introduction, the idea of maximum pseudo-likelihood is due to lBesad (|l97a l 
who first introduced the concept for Markov rando m fields in order to avoid the normalizing 
constant. This work was then widely extended and iJensen and Mailed (|l99l|) (Theorem 2.2) 
obtained a general expression for Gibbs point processes. With our notation and up to a scalar 
factor the pseudo-likelihood defined for a configuration (p and a domain of observation A is 
denoted by PLa (v?; 0) and given by 



PLa(99;0) = exp ( - / exp(-y ) TT exp (-F (xj^J \ x; 0)) 



(7) 



It is more convenient to define (and work with) the log-pseudo-likelihood function, denoted by 
LPLa (v?;6>). 



LPLa (v?; 0) = - I exp {-V {x\ip; 6)) dx - ^ V {x\^ \x;0 



(8) 



3.3 Main statistical tools 



Let us start by presenting a particular case of Campbell Theorem combined with Glotz 
Theorem that is widely used in our future proofs. For some finite configuration (resp. for some 
set G) and for all x, we denote by ^p^ (resp. Gx) the configuration (/? (resp. the set G) translated 
of x. 

Corollary 1 If the probability measure P is stationary and if the function /i(-, •) (used in Camp- 
bell Theorem) can be decomposed into h{x,ip) = l(x G A)g{x,ip) for A C IR^ where g{-,-) is 
such that g{x,ipx) = g{0,(p) for all x, then the refined Campbell theorem combined with Glotz 
Theorem allow us to obtain 



Ep[ Yl 9{x,1>\x)) = \A\Ep(^g{0,'^)eM-Vm)) 



(9) 



Let us now present a version of an ergodic theorem obtained bv lNguven and ZessinI ((19791) 
and widely used in this paper. Let D > and denote by Aq the following fixed domain 



An 



D 



< \z\ < 



D 



where for all z € IR^, \z\ = max(zi, 2:2)- 



Theorem 2 (|Ng:uven and ZessinI (|l979h ^ Let {Hg,G G Bt} be a family of random vari- 
ables, which is covariant, that for all x G IR^ , 



HgA'Px) = Hg{v), a.s. 
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and additive, that is for every disjoint Gi, G2 G Bb, 

HG1UG2 = ^Gi + Hg2, as. 

Let I be the sub— a— algebra oj T consisting of translation invariant (with probability 1) sets. 
Assume there exists a nonnegative and integrable random variable Y such that \Hg\ < Y a.s. 
for every convex G G Aq. Then, 

lim ^Hg^ = S_E(Hao\1), a.s. 

n-^+oo |G„| |Ao| 

for each regular sequence Gn U^. 



4 Consistency of the maximum pseudo-likelihood estimator 

Maximizing the pseudo-likelihood is equivalent to minimize Un{0) defined by 

We denote by 0„ = 0„ (ip) the maximum pseudo- likelihood estimator based on the configuration 
<p, alternatively defined as 

In this section, the existence of an crgodic measure is ensured, relatively to our framework, 
by Assumptions Ei, £2'°^^ and E3. The following Assumptions are needed to derive the almost 
sure convergence of this estimator. 

Ci (A„)„>i is a regular sequence of domains such that A„ IR^ as n — > +00. 
C2 For ah 6* G 0, 

V{0\-,9)GL\Pe*). 

C3 For all 6> G \ 6>* 

Pe*{W, V{0\ip;9)y^ViO\^;e*)})>0 

C4 For all 0,6' G 0, there exists c > such that P^*— almost surely, we have 

\V (0|$; e)-V (0|$; 6') \ < \\e - 6>'||'=5(0, $) (10) 

where g{-,-) is a function such that for all x, g{0,^) = g{x,^x) and such that g{0,-) G 
L\P0*). 

Remark 1 // one only assumes the existence of an ergodic measure, in particular without As- 
sumption E3 (taken into account to express Assumptions C2 and C4) then 
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» the condition C2 becomes : for all G 0, the variables V (0|-; 6) exp { — V (0|-; 0*) ) and 
exp ( — V (0|-; 9) ) are Pq* -integrable. 

• the function g{-, ■) occuring in Assumption C4 is now such that for all 6 G &, g{0,-)exp {— 
V{O\-,e))GLHP0^). 



These Assumptions have been verified in iMasa lil99d) for the Ruelle class of pairwise interaction 
function with = {(3, z) where (3 represents the inverse temperature and z the chemical potential. 

Proposition 3 Assume Pe* stationary, then under Assumptions Ci to C4, we have Pg* — almost 
surely, as n ^ +00 

e^m ^ r (11) 

Due to th e decomposition of stationary measures as a mixture of ergodic measures (see 
PrestonI (|l97a)), one only needs to prove Proposition by assuming that Pq* is ergodic. There- 
fore, in Lemmas to El Pe* is assumed to be ergodic. 

The tool used to obtain the al most sure conv ergence is a convergence theorem for minimum 
contrast estimators established by iGuvonl (|l992l ) . Define 



Kn{o,e*) = Un{o)-Un{e*) 

Lemma 4 For all 6 £ @, under Assumptions Ci and C2, we have Pg* — almost surely, as 
n +00 

Un{e) ^ U{e) = Ep^^ ( exp {-V (0|$; 6)) + V (0|$; 6) exp ( - F (0|$; 9*) )) (12) 



Proo f. Under Assumptions Ci, C2 and Ei, one can apply Theorem 01 (|Nguven and Zessin 

tozi)) to the process 

exp {—V {x\ip; 6)) dx 



H 



I An 



An 



And from Corollary^ we obtain Pg*— almost surely as n 

1 



-00 



|A„ 



Now, define 



(13) 



Let G C Aq, we clearly have 



|-?^2,g1 < ^ \V {x\^\x;G)\< E \V{x\^\x-e) 



Under Assumption C2 and from Corollary ^ we have 



^Pe* I 1^ (^1^ \x]0)\\= \MEp^^ ( \V (Ol^*; 0) 1 exp ( - y (0|$; 6**) ) ) < +00 

I 2:6* An 
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This means that for all G C Aq, there exists a random va riable Y G L^(Pf)*) s uch th at 1-^2,^1 ^ 
Y. Thus, under Assumption Ci and from Theorem [21 ( Neuven and Zessin (^3)) and from 
Corollary ^ we have P^*— almost surely 

-l-i-'i-r), A 

An 



^Ep^^ ( ^ V (xjcD \ x; 0) ) = Ep^^ (v (0|c^; 0) exp ( - F (0|c^; r) 



(14) 



We have the result by combining (|13|) and dJ. ■ 

Lemma 5 Under the conditions of Lemma^ the function Un{-) defines a constrast function, 
that is there exists a function K{-,9*) such that Pg* — almost surely the following holds for all 
6> G . 

where K(-, 9*) is a positive which, under Assumption C3 is zero if and only if 9 = 6* . 
Proof. From Lemma ^ the function K{6, 6*) > can be written 



K{0, 6*) = Epg, exp(-F (0|$; 0*)) exp(F (0|$; 9)-V (0|$; 0*))-(l+y (0|$; 0)-V (0|$; r)) 

(15) 

The result is obtained using Assumption C3 and by noting that the function t 1— > exp(f) — (1 + t) 
is positive and is zero if and only if t = 0. ■ 

Lemma 6 Under Assumption C4, the functions 9 1— > Un{9) and 9 1— > K{9,9*) are continuous 
in 9. Moreover, the modulus of continuity of Un{9) defined by 

Wn{ri) = suv{\Un{9) - Un{9')\,9,9' G0,||0-0'||<r?} 

is such that there exists a sequence {ek)k>i, with ^ as k ^ +00 such that for all k > 1 

1 



Pg* limsup Wj 



n— ++00 



0. 



Proof. Under Assumption C4, it is sufficient to prove ()16() . Denote by 
Wi,n{^) =sup||p^^ (exp(-y(xl$;0)) - exp ( - F (x|^>; 0') 



(16) 



,9,9' e &,\\9-9'\ 



and 



W- 



2,n 



i)=supJ| V {x\<^>\x;9) -V {x\<l>\x;9') 



9,9' G @,\\9-9'\\ < 



Under Assumptions E3 and C4, one can prove that P^*— almost surely 

g{x, ^)dx 



1\ expiiq 1 
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Since g{0, ■) € L^{Pg*), from Theorem 1 ( Nguven and Zessinl (1979|)) there exists Nq € IN such 
that for all n > Nq we have 

1) < '-^E,,. (,(0, *)) and W,.. (i) < '^E,,. (,(„, 



And so for all n > No, 



Wn < with 5 = AeMK)Ep^.{g{Q,^)). 



Since 

limsup I W„ ( i ] >ek\<z[^>£k\- 



Thus, it is sufficient to choose = 5'k~^ with 5' > 5 to obtain the result. ■ 
Proof of Proposition |31 

Lemma IS] and Lemma |H1 ensure the fact that we can apply Property 3.6 of iGuvonI (jlOQj) 
which asserts almost sure convergence for minimum contrast estimators. ■ 

The following proposition describes conditions C2, C3 and C4 in the case of an exponential 

functions described by ©• 



family. New conditions are denoted 03^^, Cg^^ and 04^^. For this result, let us consider energy 



Proposition 7 Conditions C2 and C4 (resp. Cs) can be replaced by (resp. C^^) where 

There exists e > such that for all i = 1, . . . ,p + 1 

Identifiability condition : There exists Ai,. . . ,^p+i, p + 1 disjoint events of Vt such 
that PQ*{Ai) > and such that for all ipi, . . . , ^Pp^i € ^1 x • • • x Ap^i the {p + 1) x (p + 1) 
matrix with entries Uj(0\<Pi) is constant and invertible. 

Proof. 

• Denote by || • ||g the norm defined for z € IR^ by \\z\\q = (X]i=i k*!'^)^'^'^ with the obvious 
notation || • || = || • ||2- We have from Holder's inequality 

|y(0|$;6>) -y (0|$;6>') | = \{9 - 0') u{0\^)\ < \\6 - 6' \\ i+e\\u{0\^)\\i+, . 

£ 

Since, is a bounded set there exists a constant k = k{£, &) such that we have 

\\e-e'\\i+.<\\e-e'\\f^ \\e - e'\\f^ < k \\e-e'\\^. 

Thus, we have (fTU)) . with c = and g{0,-) = ||u(0|<l>)||i+£. And so, 02^4" implies C4 
and obviously C2. 
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• Assumption Cg^^ means that for all y E IR^^ \ {0} and for all (/j^^, ... , (/^p+i G x • • • x 
Apj^i, the matrix (p + 1) x + 1) with entries (U)i.j = Uj (0|(/?j) (that does not depend 
on ifi, . . . , is such that U_y ^ 0. So there exists «o(2/) ^ {li • • • )P + 1} such that 

^ll^,[y\. = {p\v^o{y)) + 0" Therefore, for all y E \ {0} 

Pe^ ({(^, /w(0l99) / o}) > Pe*{A,^i^y)) > 0, 

which ends the proof. 



5 Asymptotic normality for maximum pseudo-likelihood esti- 
mates 

In this section, the existence of an ergodic measure is ensured, relatively to our framework, by 
As sumptions Ei , E!p^ and E a . The main tool used hereafter is a central limit theorem proposed 



by Jensen and Kiinsch (1994) which justifies the need of E2°'^ instead of Ei2^"^- 

To ensure the asymptotic normality for the maximum pseudo-likelihood estimator, the fol- 
lowing assumptions are needed. Denote, for some real z, by [z] the integer part of z. 

Ni The point process is observed in a domain A„ © L) = UxeA„S{x, D), where A„ C IR^ can 
de decomposed into Ujg/,jA(j) where for i = (^1,^2) 

A(i) = |z E IR2, D (^ij -^^ <zj<D (^ij - , J = 1, 2| 

for some D > 0. As n — > +cx3, we also assume that A^ IR^ such that |A„| +00 and 
\dAn\ . 

|An| " 

N2 V {0\-;9) is twice times differentiable in 9 = 6* and for all j,k = 1, ... ,p+l, there exists 
e > such that the variables 

N3 The matrix 

S(5, r ) = (^^^aI LPl£^ (^; e*)^) (17) 

is symmetric and definite positive. The vector LPL^^.'' (93; 6) is defined for any finite con- 
figuration ip and for all E and j = 1, . . . , p + 1 by 

(LPL^lj{^;0))= [ ^(x\ip;0)exp{-V{x\ip;e))dx- ^ |^(x|^\x;0). 
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N4 Vj/ G ]RP+i \ {0} 

Pe* /F«(0|v?;r)/o}) > 0, 

where for i = 1, . . . ,p + 1, (V^^HOlv, = (0|v?; r). 
N5 There exists a neighborhood V of 0* such that 1/ (•; 0) is twice times continuously dif- 
fer entiable for all j, A; = 1, . . . , p -|- 1, we have 



dV dV 



de 



< \\e-e*\\^'' /ii(o,$) 



and 



with ci,C2 > and •), /i2(-, •) two functions such that, for all x, /ij(0, $) = 
and such that /ii(0, and /i2(0, •) G ^^(Pe*). 



Jensen 



199 A) . 



Jensen and Kiinsch Hi 99. 



Re mark 2 Assump tion Ni is similar to the one of 

and lHeinricn \199^) . Among other things, Ni ensures that A„ is a regular sequence of domains 
such that A„ IE? . 



Remark 3 Similarly to Remark^ 

• the integrability condition occuring in Assumption N2 becomes : 



; e**) exp ( 



y(0|-;r))f+" and 



; e*) exp ( - y (0|-; 6*) ) are Pg* -integrable. 



• the functions hi{-,-) and /i2(', •) occuring in Assumption N5 are now such that for all 
e e V, the variables /ii(0, exp ( - F (0|-; 0) ) and /i2(0, •) exp ( - F (0|-; 0) ) are Pg*- 
integrable. Moreover, it is also assumed, for all j, k = 1, . . . ,p + 1, the Pg* -integrability of 



the variables (|^ (0|-; 9*) )^ exp ( - 1/ (0 



These Assumptions have been verified in 



Mast 



■] 0) ) and 



(Ohr)exp(-y(O!-;0) 



19991) for the Ruelle class of pairwise interaction 



function with 6 = (/?, z) where (3 represents the inverse temperature and z the chemical potential. 

For in a neighborhood of 0* , we can define, under Assumption N5, Un\9) as the vector 
derivative of Un- More precisely under Assumption Ni, we can write 

C/W(0) = |A„ri^LPLi^)(¥;;0). 



(18) 



For in a neighborhood of 0*, we can also define, under Assumption N5, the Hessian matrix 
UrS^'^ (0) given for j, = 1, . . . , p -|- 1 by 



U^n\ J A 
+ 



f dV dV 



\An\ ^ dOidOk 



89, 



{x\ip \x;6) 



_fv_ 
desOk 



{x\ip\ 0) I exp {—V {x\ip] 0)) dx 



(19) 
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Proposition 8 Assume Pq* stationary, then under Assumptions Ni to N5, we have, for any 
D fixed, the following convergence in distribution as n ^ +00 



where for some 6 and some finite configuration ip, the matrix is defined by 



%^{D,e) = \K\-'D-^Y. E 



LPL«(<^;0)LPL«((p;0) 



T 



(20) 



(21) 



'^^"\j-i\<[%]+i,j&U 



By similar arguments of lJensen and KiinschI (Il994l'l. due t o the decomposition of stationary 



measures as a mixture of ergodic measures f see iPrestonI (|l97y)), one only needs to prove Propo- 
sition IHl by assuming that Pq* is ergodic. Therefore, in Lemmas El to 1111 Pq* is assumed to be 
ergodic. 



The proof of this result is based on a general result obtained by iGuvonI ()l992l ) (Proposition 
3.7), giving conditions for which a miminum contrast estimator is asymptotically normal. The 
following Lemmas are needed to ensure these conditions. The first one ensures a central limit 
theorem for ulf'ie*). 

Lemma 9 Under Assumptions Ni, N2 and N3, 

(a) we have, for any fixed D, the following convergence in distribution as n ^ +00 



where the matrix ^{D, 6*) is defined by JiTj ). 
(6) Moreover, we have Pq* — almost surely as n ^ +00 

s„(5,r)^s(5,r). 



■p+i) 



(22) 



(23) 



Proo f, (a) The idea is to apply to Un^ {6*) a central limit theorem obtained by 

Theorem 2.1. The following conditions have to be fullfilled to apply this result. 



Jensen and Kiinsch 



i) For all i G /„ and for all j = 1, . . . ,p + 1 Ep^, ^ (LPL^^^ (^>; e*))j 



< +00. 



For all i G /„, Ep^^ ((LPL^^ («>; r)), |^.Ac^J = 0. 

\dln\ 



[iii) The set /„ is such that 



In 



0, as n 



-00. 



{iv) The matrix Varpg^ (^Kn\^/'^Un\6*)^ converges to the matrix S(D, 0), which is definite 
positive under Assumption N3. 



Condition (i) : let us write 



Ep„ 



(LPL!i)($;r)), )<2x(ri + r, 



(24) 
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where the terms Ti et T2 are respectively defined by 



f dV 

/ -T^ix\^]e*)(iM-v ix\^;e*))dx 



T2 = Epg^ 



^ |^(x|<I>\x;r 



Under Assumption N2, Holder's inequality and the stationarity of Pg*, we can prove that 



Ti < exp(3i^)|AoP x£;p^, ( [ 



dV 
90] 



dx 



< exp(3i>C)|AoP X Ep^ 
We have from Holder's inequality 



dV 



< +00. 



(25) 



T2 < Ep 



a:G<I>A 



dV 



x\^ \ x; 0*) 



And from Corollary^ it follows 

T2 < |Ao| Ep^^ ( |$aJ2 



dV 
90] 



1$;^ 



exp{-V{0\^;0*)) 



Again from Holder's inequality and under Assumption E3, we can prove that for all f] > 0, 



T2 < |Ao| eMK) Ep^^ ( I^Aol'^'^^^) Ep^^ 



dV 



3(l+»7)\ 1+^? 



Under Assumption E3 it is well-known that, for all z > 0, Epg^{\^\^^\^) < +00. Now, let e = 3ry, 
there exists k = K{e) such that 



T2 < K |Ao| Ep^ 



dV 



3+e\ l+e/3 



< +00 



(26) 



under Assumption N2. Condition (i) is obtained by combining pi]) . and 
Condition (ii) : From the stationarity of the process, it is sufficient to prove that 

EPe. ((LP<)($;0*)),|a>Ac) =0. 

Let us write for any finite configuration ip 

(LPL«(^;r)), = - / ^[x\^-e*)eM-V{x\^;e*))dx+ [ ^{x\^\x-e*)^{dx). (27) 
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Denote respectively by Gi{ip) and G2{^p) the first and the second right-hand term of ((77|) and 
by Ei = Epg^, f Gj(<I>)|<I>Ag = V'Ag)- From the definition of Gibbs point processes, 



-7—) [ <3(f^95Ao) /^2V2Ao(f^a;)lAo(x)|^(x|(^\x;0'')exp(^-y (^(^aJ(/9^c;^^ 
Denote by 99' = (v'Ao' V'Ag)- Since Q is a Poisson process we can write 
^2 = ^^/^Q(V) /^,v^'(dx)U,(x)|^(x|^\x;r)exp(-y(^,J^,.;r)) 

" Za \p^.) X ^^'^^'^ Jm:' (^I'^^o u '/'Ag \ ^; exp (-y (^:vol'^Ag; ^* 

Now, from Campbell Theorem (applied to the Poisson measure Q) 

E, = ^-J— ^ £ d. QL(V)|^ (xl^ko U ^A,; ^'^) exp {-V {^',, U xl^,.; 6^)) . 
Since from Slivnyak-Mecke Theorem, Q = Q'^, one can obtain 

= Q{d^') d. 1^ {.y^^ u ^As; exp (-y {^'^^ u xI^a^; ^'^)) 

= — r/ QidifAj [ dx^ix\ip;e*)exp{-V{xy,e*))exp(-V (ipj^Jip^c;e 



Condition {in) : this condition is equivalent to Assumption N 

(1) 



Condition {iv) : let us start by noting that the vector LPlI^'' {ip; 0*) depends only on 99^ for 



j such that \j — i\ < 



+ 1. Prom (jlHp . we can obtain 



Varg. (|A„ti/2[/(i)(0-)) = |A„|-iVarp^. (LPLJi^ ($5^) 



lA„ri ^ £;p,. (LPLi^)(c^;r)LPL^)(cI>;rn 

+ J] £;p,. (lplI^) (CD; r ) LPL« (CD; r f ) I . 



Let j € /„ such that [j — 'i| > 



+ 1, then using condition (ii) 



Ep^, [ LPl!^^) (CD; r) LPL^') (cD; 0*)] = Ep^, (e ( LPlJ^^ (cD; 0') LPL^'^ (cD; e*f\LFL^l^, (cD; 6*) 



Ep^, (e (LPLi^) (CD; ^ ) |LPL« (cD; 0*)) LPL« (cD; 0^) 
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Now, denote by / the following set 

I = < k e In,\k - i\ < 



+ l,Vi e5/„ 



and (for the sake of simplicity) by Eij the following mean 



(1) 



From the stationarity of the process, we can write 



Varp,. (|A„|V2t7(i)(r)) = \A^\-' 



\ 



E E ^M + E E 

yiGln\I |j_i|<[£j+l,jG/„ \j^i\<[%]+l,jeln J 

= ii„\7| |A„|-i Yl ^0.^ + 1^1 i^nr' E ^^oj> 

N<[§]+i b--ol<[§]+i 
for some iq G din- From the definition of the set /„, we have as n — > +00 

Varp,. (|A„r/2[/(i)(0^)) ^ ^ Eo, = ^[5,0*) . 

Kl<[§]+i 

(6) According to (PT|). it is easy to see that ^^(D, 0*) is defined such that as n ^ +00, 



We lea ve the reader to check that under Assumption Ni and from Theorem 1 (|Nguven and Zessin 
(19791)), we have P^^-almost surely as n ^ +00, S„(5, 6*) 6*). ■ 



Remark 4 From the previous proof, we can note that Assumption N3 is fuUfilled as soon as one 
can prove that for n sufficiently large the matrix Yarp^^ MA„|^/^?7n^^(0*) ] is definite positive. 



Lemma 10 Under Assumptions Ni, N2 and N5, there exists a neighborhood V of 9* on which 
Un{-) is twice times continuously differentiahle and a random variable Y such that for all j, k = 
1, . . . ,p + 1 and for all 9 we have, 



uJ^H9- 



< Y. 



Proof. Let j,k = 1, . . . ,p-\-l. Under Assmuption N5, there exists a neighborhood V of 9* such 
that we can write for any configuration ip 

1 f d^V 



Uj'\9))^ 



\^n\ Ja„ dOjdOk 
dV 



{x\(p] 9) exp {—V 9)) dx 
dV 



1 f dV dV 

I n\ J j k 



1 ^ 



{x\ip \ x;9) . 



(28) 
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Denote respectively by Ri,R2,Rz the three right-hand terms of the previous equation. Under 
Assumption N5, one can choose the neighborhood V such that \\6 — 6*\\ < k. Thus, one can 
obtain 



h2{x,ip) + 



dOjdek 



|fllI<e>:pW^/^_(.' 

\R2\ < exp(K)-^ / (K^^'hi{x,ipf + K''hiix,^) 



dV 



dx 



:x\^;e*) 



+ K^^hiix, If) 



dV 



ix\ip;e'') 



+ 



l^sl^hO ^ U'''h2ix,^\x) + 



dV dV 
-(.|.;«*)_(.|^;«.) 



Under Assumptions Ni and N2, from Theorem [21 (|Nguven and ZessinI ()l97fll )). and using the 
stationarity of P^*, there exists Nq G IN such that for all n > Nq, we have P^*— almost surely 



dOide 



•jOVk 



|Pl| < 2 X Q-K.Y>{K)Ep^, (^K^2/j2(o,$) + 
|i?2| <2xexp(ir)|£;p,. |^fe2-i/,,(O,^.)2 + ||^(O|$;0 



dV 

del. 



dV 



d9i 



IP3I < 2 X exp{K)Ep^, K^2/i2(0,^>) + 



Consequenlty, for n large enough, there exists a positive constant k' such that y Un ^'^^ {0) ) 
n' , which implies the result. ■ 

Lemma 11 Under Assumptions Ni and N2, we have almost surely, as n ^ +00 
where t/^^^(0*) is the (p + l) x {p+ 1) matrix whose entry 



j,k 



< 



IS 



dV 



dV 



|«>;r)exp(-y(0[$;r)) 



Furthermore, under Assumption N4, U_^'^^ is a symmetric definite positive matrix. 



(29) 



Proo f. Let j,k = 1, . . . ,p+l. Under Assumptions Ni and N2 and from Theorem01 (|Nguven and Zessin 
we have almost surely, as n ^ +00 



|Ac 



e") exp {-V (x|$; 6*)) dx 



+ ^^Ep^* 1^ {A^-^ 0*) ^ (xjcK; r) exp {-V (x|<I>; r)) dx^ 



jdOk 



{x\^\x;e*) 



(30) 
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Equation (|29() is obtained using Corollary ^ And under Assumption N4, it is easy to see that 
is a symmetric definite positive matrix. ■ 



Proof of Proposition IS] Using Lemmas El a 1111 one can apply a classical re sult concerniii g 
asymptotic normality for minimum contrast estimators, e.g. Proposition 3.7 de iGuvonI (|1992t ). 
in order to prove as n ^ +00 

|A„|V2s„(Ar)-i/2[/(2)(r) - a/- (0,4+1)- 



The result is then obtained using the fact that 0„ is a consistent estimator of 0*. ■ 

Let us precise, as in Section ^ the different Assumptions for energy functions that can be 
written as @- 

Proposition 12 For energy functions described by Assumptions N2 and N5 (resp. 'N4) 
can be replaced by N2''^ (resp. N^'^^j 

^2^5^ ^^''^ i = 1, . . . ,p + 1, there exists e > such that nj(0|-) G L^"*""^ (Pg* ) . 

-I ex 
'3 



^exp ^ (^exp 



The proof is trivial. 



6 Some examples 

In this section, it is assumed that the sequence of domains satisfies Ni (which implies Ci). 
Moreover, we only focus on examples satisfying the following convenient Assumption denoted 
by M : 

M There exists Ki,K2 > such that for any finite configuration ip, we have for all x 

-Ki < Ui{x\ip) < K2, for z = 1, . . . ,p + 1. 

Quite obviously. Assumption M ensures 02^4" and N2''^. Let us now present a Gorollary of 
Propositions |21 and |H1 

Corollary 13 Under Assumption M and , the consistency of the maximum pseudo-likelihood, 
that is the result Ul\) . is valid. And in addition with Ng''^, its asymptotic normality property, 
that is the result \2U\) . is ensured. 
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6.1 Pairwise /?-Delaunay model 



We first deal with our main example : /3-Delaunay (of order some small enough fixed Pq) 
model with multi-Strauss pairwise interaction function. In other words, 



(31) 



with ui{ip) = \ip\ and for any i £ {2, . . . ,p + 1}, 



Ui{ip) 



E 



'-]di-i,di] 



where = di < d2 < ■ ■ ■ < dp+i are some fixed real numbers. Literally, Ui{ip) {i > 1) corresponds 
to the number of (/3-Delaunay) edges of length between and di. We may also notice that the 
range of the pairwi s e inter action function is dp+i, that is ii(^)(^; 93, 0^^^) = when ||^|| > dp^i. 

In iBertin et al.l (|l999a|), it is proved that this model satisfies Assumption M. Let us now 
verify the technical conditions Cg^^ and Ng^^. 



Proposition 14 Assumption is satisfied for the (3-Delaunay model with multi-Strauss 

pairwise interaction function. 

Proof. Denote by A the following domain 

A = {zeM^ ■.-D<Zi <D,i = 1,2} 

and by Ai the event Ai = {ip, ip^^ = 0}. We clearly have for all S Ai, u{0\ipi) = (1,0,..., 0)'^. 
Now, let us give for j = 2, . . . ,p + 1, the points cij and 02 j such that the distances d{0, cij) = 
d{0, C2j) = d{cij,C2j) = ^^i^it^i. Denote for j = 2, . . . ,p + 1 the following events for some 77 > 

Ajiv) = \^(p & ^ ■■ (PA = {zi,Z2},zi £ B{cij,7]),Z2,j G B{c2,j,r])y 

One can choose rj such that for all if G ^j(??), the distances d{0, zi),d{0, Z2 and d{zi,Z2) are 
comprised between dj-i and dj. One can also choose rj such that the smallest angle of the 
triangle with vertices {0,cij,C2j} is strictly greater thant /?q, which means that {0,cij,C2j} € 
Del^°{<Pj). Now, it is easy to see that the matrix U_ defined in Proposition [3 is given by 



U = {uj{0\ipi))i<ij<p+i 



/l 

1 3 

1 

1 



\ 




3 



and is clearly invertible, which ends the proof. 
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Proposition 15 Assumption N3 is satisfied for the /S-Delaunay model with multi-Strauss pair- 
wise interaction function. 

Proof. From Remark0J it is sufficient to prove tliat the matrix Yarp^^, {\An\^/'^Un\^', 0*)) is 
definite positive f or n sufficientlv larg e . Let D > D, y ^ IR^^^ and let A = U|j|<iA(j), by the 
same argument of l.Tensen and KiinschI ()1994l ) (Equation (3.2)), we can write 

/Varp^. (|A„|i/2[/(i)) y > |A„! Ep^^ (varp,. (/[/^ |cI>a^, £ ^ 37Z^)) . 

Now, following the proof of Lemma|Hl((a) condition (iv)), one can prove that there exists ng G IN 
such that for all n > uq, 

/Varp,. (|A„[i/2) y > ^Ep^^ (varp,. (/LPL^^^ ($;0*) |$a,, 1 < \£\ < 2)) . 

The aim is to prove that the function h{^) = y^LPL^3^ ($; 0*) is not almost surely a constant, 
when the variables ^a^, 1 < K| < 2 are (for example) fixed to 0. Assume that the function h(-) 
explicitly given for any finite configuration ip by 

K'P) = '^ykl L Ukix\iPAjexp (^-0*^u(x|(/9ao)) dx - Y1 ^fe(^l<^Ao \x) 
k=i [-^^ xe^^^ 

is constant for ah 99 G fi' = {(/? e : ipj^^ = 0, 1 < |^| < 2}. 

Denote by Aq = {99 G il' : c^Aq — ^} ^^'^ W ^1 — {9? G Jl' : Iva,,! — clear that, 

P0*{Ao) > and P0*{Ai) > 0. We have for all (pQ G Aq and for all (pi G Ai 

hi^o) =yi\A\exp{-el) and h{ip^) = yi\A\exp{-9{) - yi. 

Assuming h(-) constant implies that yi = and then h{-) vanishes. We now consider particular 
configurations of two points in Aq and empty in A^ \ Aq, 1 < |^| < 2. Let us first introduce the 
following sets for any j G {1, ■ ■ ■ ,p — 1} and any rj > 



D,{r,) = {(zi,Z2) G A2 : zi G e ((0,0), |) and Z2 ^ B (^(d„0), } 
D-{i^) = {(zi,Z2) G A2 : G e ((0,0), |) and e B {{d, - |,0), |) } C D,(r?) 
Z)+(7?) = {(zi,Z2) G A2 : zi G e ((0,0), |) and ^ B ((d, + |,0), |) } C D,(r?) 

When 7/ is small enough, the couple of points (2:1,22) G Dj{r]) (resp. Dj{r]) and D^{r])) are 
such that < dj — r] < d{zi, Z2) < dj + r] < dj+i (resp. < dj — r] < d{zi, Z2) < dj and 
dj < d(zi, ^2) < dj + rj < dj+i). 

We now derive the corresponding events for any j G {1, ■ ■ ■ ,p — 1} and any > 

Ajiv) = jy? G : (/?Ao = {^1,2:2} with (21,22) G Dj{r])^ 
Ajiv) = {^en' : (^Ao = {^1,^2} with (21,22) G DT{ri)^ c ^j(r?) 
^/(??) = G : (^Ao = {^1,^2} with (21,22) G D+{7])^ C Ajirj) 
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satisfying Pe*iAj{r])) > 0, Pe*iAj{7])) > and P0*iA+{r])) > 0. 

Let us fix some ip € Aj{r]). There exists some unique couple of points (21,^2) G Dj{''l)i for 
which we define the fohowing domain 

A(zi,2;2) = jx e A : {x,zi,z2 }GDe/f°(^U{x})}. 

Since {zi,Z2} ^ Del2°{^p)-, we then derive that 

= h{(p) = ^ykl Uk{x\{zi,Z2})exp(-0*'^u{x\{zi,Z2})jdx. (32) 



When {x, zi, Z2} G Z?e/3°((/?U{3;}), we decompose Uk{x\{zi, Z2}) into two additive terms in order 
to isolate the contribution of x : 

Uk{x\{zi,Z2}) = Z2) + ul{zi,Z2) 

with u^^{zi,Z2) = and uf{zi,Z2) = 1, and for k ^ 1, u~''{zi,Z2) = - Z2II) and 

-"^(21,^2)= Yl hdk-iMi\\x- ZjW)- 

Then Equation ()32|) becomes for any y € IR^^^ with yi = 



h{(f) = = exp\^-e*^u ''{zi,Z2)j {ifu ^(zi, 2:2)/! (21,^2) + 2/^/(21,^2)) 
^ /w-"(2l, Z2)flizi,Z2) + y^f{zi,Z2) = (33) 

where u"^ = (uf^, • • • , u~^-^), tx^ = (nf , ■ • • , u^.^^) and / = (/i, /2, • • • , /p+i) with 

fk{zi,Z2) = I ul{zi,Z2)exp(-0*'^u''{zi,Z2))dx 

Since each /fc is continuous, one could assert that for any e > 0, there exists 7/ > such that 
for any {z\,Z2) G Djij]), \fkizi,Z2) - fk\ < e where fk = /fc((0, 0), (dj, 0)) is positive. We then 
set 5k{zi, Z2) = fk{zi, Z2) — fk and S = {5i, - ■ ■ , (5p+i). We now apply the equation (|33|) for some 
fixed ifj G Aj{r]) and 93+ G ^j'^(^)- By denoting (zf ,^2^) G Dj{'q) and (zj'',^:^) G Dfit]) such 
that n Aq = {zjf, ^2"} and n Aq = {z^ , z^}, we have 

yj/i(2r'^2") +'!/'f{Zl,Z2) = 
y,+i/i(z+, Z2+) + y^fizf, z+) = 

By substracting these two terms, we can obtain 

{Vj+i - yj)h = yjh{zi,z~) - yj+i(5i(2;+,z+) {6{z~ ,z~) - 5{z^,z^)) (34) 



By the previous continuity argument on the fk, on can choose > (depending on y) small 
enough such that the absolute value of the right-hand term of (|34|) could be lower than any 
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e > 0. Thus, by assuming that yj / Vj+i and choosing e = ^\yj+i — Ujlfi, there exists rj such 
that \yj^i—yj\fi < ^\yj^i — yj\fi which leads to an obvious contradiction. Thus, (|SH) holds only 
if yj+i = yj- By iterating this argument, we obtain that y2 = Vs = ■ ■ ■ = Up+i and by applying 
this result on the equation (|33() . one may assert that for any y € IR 

P+i 

y^fk{zi,z2) = 0, 

k=l 

which implies that y = since fk{zi, Z2) is positive for any {zi, Z2) € Dj{ri) {j arbitrarily chosen 
in{l,---,p+l}) ■ 

We propose a simulation study to verify the consistency of maximum pseudo-likelihood es- 
timator. We consider the model (|31() with the vector of parameters 9 = (0,2,4). The vector 
of bounds d is assumed to be known and fixed to d = (0,20,80). The simulation proce- 
dure used here is a dire c t ada p tation t o the Delaunay energies of th e Geyer and M0ller pro- 



posal (|Gever and Moiled ((1994 ). iGeved (jl999D). We refer the reader to iBertin et al.l (|l999dh for 
a detail of the used algorithm. One simulation of such a point process is proposed in Figure ^ 
Table n summarizes the different results obtained via m = 5000 replications each one is gener- 
ated after one million of iterations of the algorithm. One may verify that both the bias and the 
standard deviation become smaller and smaller as the domain A.„ grows. 



6.2 Other examples of pairwise interaction models 

In order to satisfy Cg^^ and Ng^^ for models on the complete graph or on th e k nearest 



Jensen and Kilnsch 



neighb ours graph with multi-Strauss pairwise interaction function, we can chosse as in 
( 1994} ) a configuration with one point or two points. On the delaunay graph, it may be interest- 
ing to study multi-Strauss interaction function on the circumradius or on the smallest angle of 
each Delaunay triangle. As discussed previously the identifiability assumption holds easily 
but Ng''^ needs more attention. Otherwise, for pairwise Delaunay model, we can replace the as- 
sumption on the smallest angle by a hard-core assumption and then Q, by the set of admissible 
configurations = {99 G ri : Vx, y€ipxip,x^y \\x — y\ \ > 5}. 
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Fig. 1 - Example of the points, with (left) and without the edges (right), of the realization of 
the /3-Delaunay model with multi-Strauss pairwise interaction function where parameters 6 and 
d are respectively fixed to (0, 2, 4) and (0, 20, 80). 



Domain A„ 


Mean of Estim. of 6*2 (Std Dev.) 


Mean of Estim. of 6*3 (Std Dev.) 


[-250,250]2 
[-350,350]2 
[-450,450]2 


2.068 0.104 
2.049 0.071 
2.041 0.056 


4.382 0.786 
4.223 0.551 
4.144 0.436 



Tab. 1 - Empirical mean and standard deviation of maximum pseudo-likelihood estimates of 
parameters of 02 = 2 and ^3 = 4 representing the levels of a multi-Strauss pairwise interaction 
function where the vector of bounds is assumed to be known and fixed to d = (0, 20, 80) . These 
results are obtained from m = 5000 replications of the point process described by ()31() generated 
in the domain [—600,600]^. Three sizes of domains A„ have been considered. 



